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ABSTRACT 


During  the  course  of  this  program  experimental  investigations  have 
been  carried  out  to  meet  two  objectives,  a)  to  measure  the  temperature 
dependence  of  rate  coefficients  of  radical -molecule  and  molecule-molecule 
reactions  of  specific  importance  in  modeling  of  combustion  and  plume  processes, 
and  b)  to  contribute  to  the  general  understanding  of  the  temperature 
dependence  of  bimolecular  reactions,  particularly  as  it  relates  to  observed 
curvature  in  Arrhenius  plots.  The  technique  of  flash  photolysis-resonance 
fluorescence  has  been  adapted  to  measure  directly  rate  coefficients  for  OH 
reactions  in  the  temperature  range  of  250-1100  K.  Kinetic  data  for  the 
important  reactions  of  OH  +  Mg,  Dg,  CH^,  CgHg,  CgHg,  HC1,  and  CO  was  obtained. 
In  addition,  the  general  oxidative  conversion  of  BClg  and  BFg  were  also 
investigated. 


INTRODUCTION 


Aircraft  and  rocket  exhaust  plumes  provide  an  energetic,  highly  non- 
uniform,  local  environment  within  which  microscopic  physical  and  chemical 
processes  act  to  dictate  the  nature  of  bulk  observations  of  the  medium. 
Analyses  of  fuel  performance  and  tracking  procedures  rely  on  modeling  of 
those  contributing  processes  which,  in  combination,  govern  the  efficiency 
of  the  propulsion  system  or  surveillance  method  under  examination.  Chem¬ 
ically,  reactions  within  the  combustor  and  afterburning  within  the  plume 
place  heavy  demands  on  the  estimating  prowess  of  kinetic  modellers.  In 
most  cases,  it  is  simply  not  possible  to  calculate  which  among  the  myriad 
of  possible  processes  are  most  important  (and  thus  which  among  the  possible 
molecular  species  are  of  greatest  concern)— an  extensive  data  base,  promot¬ 
ing  realistic  modeling  parameter  limitations,  is  needed.  From  a  kinetic 
standpoint,  optimized  acquisition  of  this  data  base  implies  the  measure¬ 
ment  of  oxidant/fuel/reaction  intermediate  reaction  rate  constants  needed 
directly  for  model  input  and/or  required  as  bases  for  the  generation  of 
'global'  kinetic  equations  applicable  to  entire  families  of  chemical  re¬ 
actions.  Improvement  of  the  understanding  of  those  factors  which  govern 
the  temperature  dependence  of  combustor/plume  relevant  bimolecular  reac¬ 
tions  was  the  primary  goal  of  this  research  program. 

The  results  are  presented  in  the  following  four  Sections.  Section  I 
has  been  recently  published  in  the  Journal  of  Physical  Chemistry.  Section 
III  will  be  submitted  shortly  to  the  Journal  of  Physical  Chemistry  for 
publication.  It  Is  anticipated  the  contents  of  Section  IV  would  also  be 
submitted  for  publication  in  the  next  few  months. 


I.  A  FLASH  PHOTOLYSIS-RESONANCE  FLUORESCENCE  KINETIC 
STUDY  OF  THE  REACTIONS  OH  +  H£  -  HgO  +  H  AND 
OH  +  CH4  H20  +  CH3  FROM  298-1020  K 


(Published  in  the  Journal  of  Physical 
Chemistry,  8£,  3126  (1980)) 
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INTRODUCTION 

Experimental  measurements  of  absolute  rate  constants  for  radical - 
molecule  reactions  have  historically  been  divided  according  to  tempera¬ 
ture  regime.  At  low  temperature  (T  <  500K)  a  variety  of  techniques*  have 
.“been  developed  and  utilized  in  recent  years  in  response  to  the  input  demands 
of  atmospheric  chemical  modeling.  Generally  speaking  these  techniques  per¬ 
mit  the  adjustment  of  experimental  conditions  such  that  the  absolute  rate 
constant  measurements  for  many  specific  radical-molecule  reactions  may 
be  made  with  little,  and/or  quantifiable,  interference  from  competing  re¬ 
action  processes.  At  high  temperature  (T  ;>  1000K)  kinetic  isolation  of 
a  given  reaction  has  proven  to  be  much  more  difficult  to  achieve.  Experi¬ 
ments  at  flame  temperatures  typically  range  from  studies  utilizing  end 
product  analysis  alone  to  investigations  involving  detailed  concentration 
mapping  vs.  time  (distance)  of  a  number  of  radical  (stable  molecule) 
species  whose  interactive  chemistry  is  very  complex.  Sequences  of  re¬ 
action  schemes  with  parametrized  absolute  rate  constant  values  are  used 
to  iteratively  reproduce  the  measured  concentration  profiles  (product 
yields).  Frequently  the  extracted  rate  constants  are  model-dependent, 
ratioed  to  those  of  another  insufficiently  characterized  reaction,  or 
sensitive  within  only  broad  limits  on  the  obtained  raw  data.  The  inter¬ 
pretive  clarity  of  high  temperature  kinetic  measurements  has  been  observed 
to  parallel  the  level  of  experimental  control. 

At  temperatures  intermediate  between  these  separated  regimes,  500- 

1000K,  relatively  few  investigations  of  radical -molecule  reaction  rate 

constants  have  been  undertaken.  Two  reactions  which  have  received  some 
2  3 

attention  *  in  this  interval  as  well  as  extensive  study  at  high  and 
low  temperatures  are  the  subject  of  the  present  investigation: 
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OH  +  H2 


H20  +  H 


(1) 


k2 

OH  *  CH„  -A-  H20  +  CHj  (2)  # 

Published  kinetic  data  on  Reaction  (1)  has  recently  been  critically 
m  -  4 

-  evaluated  and  summarized  by  Cohen  and  Westberg.  Based  principally  on 
kinetic  measurements  made  below  500K  and  above  1000K,  these  authors 
recommend  a  three-parameter  rate  constant  expression  of  functional  form 
k(T)  =  ATnexp-(E0/RT)  as  an  appropriate  fit  to  the  Arrhenius  graph  cur¬ 
vature  experimentally  defined  for  Reaction  (1).  Non-Arrhenius  behavior 

c 

is  also  evident  in  kinetic  measurements  made  on  Reaction  (2).  Zellner 
has  used  this  same  functional  form  to  empirically  effect  a  smooth  joining 
of  the  low  and  high  temperature  rate  constant  measurements  for  this  reaction. 

In  this  paper  we  present  the  first  in  a  series  of  studies  of  the 
reaction  kinetics  of  the  hydroxyl  radical  at  intermediate  temperatures. 

Using  the  direct  kinetic  technique  of  flash  photolysis-resonance  fluor¬ 
escence  absolute  rate  constants  for  Reactions  (1)  and  (2)  have  been  ob¬ 
tained  up  to  temperatures  of  1000K.  These  results  will  be  described  and 
compared  with  those  of  previous  measurements  and  compilations. 

EXPERIMENTAL  SECTION 

The  application  of  the  flash  photolysis-resonance  fluorescence  tech¬ 
nique  to  the  study  of  OH  radical  reaction  kinetics  has  been  discussed 
previously.^  Hence  we  will  describe  the  experimental  apparatus  briefly, 
emphasizing  only  those  new  features  which  for  the  first  time  extend  the 
temperature  range  of  applicability  of  the  technique  to  above  1000K. 

A  schematic  diagram  of  the  experimental  apparatus  is  shown  in 
Figure  1.  The  principal  components  of  this  system  are  (1)  a  quartz 
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reaction  cell  resistively  heated  using  electrically  insulated  tantalum 
wire  windings  mounted  to  its  graphite-coated  outer  surface;  (2)  a  spark 
discharge  flashlamp  perpendicular  to  one  face  of  the  cell;  (3)  a  CW  OH 
resonance  lamp  perpendicular  to  the  photolysis  beam;  (4)  a  photomultiplier/ 
bandpass  filter  combination  for  monitoring  OH  resonance  fluorescence 
perpendicular  to  both  the  photolysis  and  resonance  radiation  beams;  and 
(5)  a  signal  averager  and  fast  photon  counting  electronics. 

The  present  experiments  were  carried  out  using  a  static  reactor  con¬ 
figuration.  Typical  gas  mixtures  consisting  of  150  mtorr  H20,  0-1  torr 

MCH„),  and  50  Torr  argon  were  flash  photolyzed  thereby  initiating  the  pri¬ 
mary  reaction  sequence 

H2o - ^ - >  H  +  0H(X2ti) 

165  nm  <  X  < 185  nm 

kM 

OH  +  H2(CH„)  -£!  H20  +  H(CHj). 

It  is  worth  noting  here  that  our  use  of  Suprasil  quartz  optics  with 
a  transmission  cutoff  at  165  nm  prevented  photolysis  of  H2  and  CH.,^  and 
thus  eliminated  a  major  potential  source  of  unwanted  radicals  and  com¬ 
peting  radical-radical  reaction  processes.  Following  the  flash  weakly 
focused  resonance  lamp  radiation  continuously  excited  a  small  fraction 
of  the  OH  to  the  A2Z+  state  and  the  resultant  (0,0)  band  fluorescence 
was  counted  in  real  time;  the  temporal  profile  of  OH  decay  was  constructed 
by  signal  averaging  over  a  large  number  of  flashes  for  each  fixed  reactant 
concentration.  The  reactant  concentrations  [H2]  and  (CHh ]  were  maintained 
at  large  excess  over  l  OH]  (lOHl^.g  =  1-5  x  1011  molec/cm* )  and  the  reaction 
kinetics  were  thus  pseudo  first  order  in  10H] .  Exponential  [OHj* 
decays  (through  several  1/e  times)  were  obtained,  their  time  constants 
yielding  the  first  order  rate  constant  for  OH  disappearance  at  the  re- 
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actant  concentration  of  that  mixture.  Bimolecular  rate  constants  were 
obtained  as  the  slope  of  a  plot  of  the  measured  first  order  rate  constant 
vs.  reactant  concentration  ([H2]  or  t CHW 3 ) .  System  parameters  such  as 
flash  energy,  H20  concentration  and  total  pressure  were  varied  to  demon- 
-“"strate  the  lack  of  importance  of  competing  secondary  processes  in  our 
measurements  on  Reactions  (1)  and  (2).  During  the  generation  of  most  10H] 
decay  curves  several  cycles  of  reactor  filling/evacuation  with  identical 
gas  mixtures  were  used  so  as  to  negate  the  kinetically  deleterious  effects 
of  reactant  depletion  and  product  buildup  (within  the  exposure  limits 
used  the  measured  rate  constants  were  shown  to  be  independent  of  the 
number  of  flashes  to  which  any  gas  filling  was  subjected). 

The  two  principal  experimental  factors  which  have  permitted  direct 
measurements  of  absolute  rate  constants  to  be  made  at  temperatures  in 
excess  of  1000K  involve  the  system  optics  configuration  and  the  reaction 
cell  design.  Collimation  and  weak  focusing  of  the  incoming  flashlamp 
and  resonance  lamp  radiation  and  of  the  emitted  fluorescence  have  been 
employed  and  have  resulted  in  an  OH  detection  sensitivity  which  is 
significantly  enhanced  relative  to  that  available  in  previous  high-pressure 
studies  of  Reactions  (1)  and  (2).  A  detailed  discussion  of  the  detection 
sensitivity  characteristics  of  a  completely  analogous  experimental  con- 

O 

figuration  is  contained  in  a  recent  publication  describing  our 

laboratory's  study  of  the  OH  +  NO*  +  M  HNO3  +  M  reaction.  Highly 

polished  Suprasil  windows  have  been  fused  using  quartz  tape  to  an  all 

g 

quartz  body  to  form  the  high-temperature  reactor.  Efficient  heating 
and  radiation  shielding  maintain  a  higher  temperature  at  the  body 
of  the  reactor  than  that  at  the  protruding  windows,  the  reaction  tempera¬ 
ture  being  measured  by  both  internal  (insertable)  and  external  (body- 
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mounted)  Chrome! -Alumel  thermocouples.  During  the  period  of  measurement 
of  each  individual  bimolecular  rate  constant  the  reactor  temperature  was 
held  constant  to  within  ±  3K  by  applying  a  controlled  constant  voltage  to 
the  resistive  heaters.  The  temperature  profile  along  the  axis  of  the 
cylindrical  reactor  was  measured  using  the  retractable  thermocouple  at 
pressure  conditions  identical  to  those  existing  during  kinetic  measure¬ 
ments.  The  temperature  differential  between  the  upper  and  lower  portions 
of  the  reaction  zone  (maximum  length  2  cm)  was  found  to  be  *  5K  at  a 
reactor  temperature  of  lOOOK.  The  bottom  window  (fluorescence  port) 
temperature  was  measured  to  be  %  25K  below  that  of  the  reaction  zone  at 
1000K;  given  the  heater  and  radiation  shield  positions  and  the  overall 
reactor  geometry,  side  window  temperatures  are  estimated  to  be  comparable 
to  that  measured  at  the  bottom  window.  Considering  the  relatively  large 
reactor  diameter  used  in  this  work  (reactor  volume  %  500  cm3),  such 
temperature  gradients  are  likely  to  be  only  a  minimal  source  of  experi¬ 
mental  error. 

Figure  2  shows  typical  OH  decay  profiles  at  three  temperatures. 

These  decay  rates  were  all  measured  using  the  same  diluent  gas  pressure, 

H20  pressure,  and  photolysis  flux.  The  high  temperature  experiments 
required  longer  periods  of  signal  averaging  than  the  lower  temperature 
measurements;  for  example,  twice  as  many  flashes  were  needed  at  1020K 
than  at  619K  to  produce  a  good  quality  decay  curve.  At  an  H20  partial 
pressure  of  150  mTorr  the  fluorescence  signal  strength  as  a  function 
of  temperature  was  seen  to  scale  closely  with  H20  number  density  indicating 
minimal  temperature  dependent  optics  transmission  effects;  repeated 
temperature  cycling  of  the  reactor  did,  however,  induce  long-term  degrada¬ 
tion  of  the  vacuum-UV  transmission  efficiency  of  the  Suprasil  windows. 
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[The  static  reactor  used  in  the  present  work  has  since  been  replaced  by 
an  analogous  "slow  flow"  high- temperature  quartz  reaction  cell  which 
transmits  a  continuously  controlled  flow  of  mixture  of  desired  composition 
and  thereby  permits  each  initiating  flash  to  encounter  a  "locally  fresh" 
reaction  mixture.] 

The  gases  used  in  this  study  were  obtained  from  Matheson  Co.  and 
had  stated  minimum  purities  as  follows:  H2,  99.9999*;  CH.,,  99.99*; 
and  Ar,  99.9995*.  Prior  to  their  use,  hydrogen  and  methane  were  passed 
through  a  trap  at  77K  to  remove  residual  condensibles. 

RESULTS  AND  DISCUSSION 

The  bimolecular  reaction  rate  constants  obtained  in  this  study  were 
generated  from  measurements  of  hydroxyl  radical  decay  rates  at  several 
fixed  reactant  concentrations  for  each  temperature.  Plots  of  the  first 
order  decay  constant  vs.  reactant  concentration  at  selected  temperatures 
between  499-1020  K  are  given  for  Reactions  (1)  and  (2)  in  Figures  3  and 
4,  respectively.  The  bimolecular  rate  constants  obtained  as  the  slopes 
of  the  included  least-square  lines  are  listed  with  ±2o  precision  limits 
in  Tables  I  and  II.  Room  temperature  rate  constant  values  obtained  in 
this  work  are  in  good  agreement  with  those  of  previous  studies  tabulated 
in  the  data  evaluations  of  References  4  and  5. 

The  major  potential  sources  of  systematic  error  in  the  measured 
values  of  kj  and  k2  are  the  uncertainties  in  the  reactant  concentrations 
and  in  the  temperature.  At  298K  the  reactant  concentrations  of  source 
mixtures  could  be  determined  quite  accurately  (%  3%)  using  a  capacitance 
manometer  since  both  H2  and  CHi,  are  non-sticky  and  non-condensible  com¬ 
pounds.  Reaction  fillings  were  made  using  these  calibrated  source  mixtures 
and  elevated  temperature  reactant  concentrations  were  calculated  by  scal- 
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ing  accordingly  to  source  mixture  content  and  temperature.  Including 
errors  due  to  temperature  uncertainties,  then,  the  estimated  error  in 
reactant  concentrations  was  *  10*.  Thus  we  estimate  that  the  accuracies 
(95%  confidence  limit)  of  the  measured  values  of  ki  and  k2  are  n,  13% 
at  29SK  increasing  to  20%  at  1000K. 

Our  data  for  ki(T)  are  plotted  along  with  selected  previous  results 
and  analytical  evaluations  in  Arrhenius  graph  form,  £n  k(T)  vs.  1000/T, 
in  Figure  5.  Our  results  cover  a  sufficiently  wide  temperature  range  to 
independently  demonstrate  the  existence  of  curvature  in  the  Arrhenius 
plot  for  Reaction  (1).  At  low  and  relatively  high  temperature  our  measure¬ 
ments  are  in  excellent  agreement  with  the  rate  constant  values  calculated 
from  Cohen  and  Westberg’s  summary  expression^  ki(T)  =  1.83  x  10-ls  T13 
exp(-1835/T)  cm3molec“1sec_1.  In  the  interval  499-648  K,  however,  measure¬ 
ments  made  in  the  present  study  fall  about  20%  below  those  derived  from 
this  ki(T)  expression.  Similarly,  in  the  only  previously  published  direct 

kinetic  study  of  Reaction  (1)  at  intermediate  temperatures,  Westenberg  and 
2 

deHaas  measured  ki  values  between  403-745  K  which  averagd  22%  below  those 
calculated  from  the  Cohen  and  Westberg  formula.  Interestingly  Cohen  and 
Westberg  did  not  weight  data  at  intermediate  temperatures  heavily  in  their 
evaluation;  indeed,  hardly  any  measurements  between  500-1000  K  were  avail¬ 
able.  While  our  measurements  yield  generally  very  good  agreement  with 
the  ki(T)  values  calculated  from  Cohen  and  Westberg's  expression,  they  do 
suggest  somewhat  more  pronounced  Arrhenius  graph  curvature  than  that 
expression  describes.  The  best  fit  expression  calculated  from  our  data 
only  is  ki(T)  =  4.12  x  10"19  T2**  exp(-1281/T)  cm3molec"1sec”1 .  Zellner's 
suggested^  expression  of  ki(T)  *  1.66  x  10”ls  T1,6  exp(-1660/T)  cn^moleiT'sec 
is  perhaps  the  most  appropriate  fit  to  all  available  data.  [With  regard 
to  curvature  in  Arrhenius  plots  it  is  of  interest  to  note  that  preliminary 


results  of  our  OH  +  02  -*■  DHO  +  D  study^  between  298-932  K  demonstrate 
markedly  less  Arrhenius  graph  curvature  than  that  found  to  be  appropriate 
for  Reaction  (1).] 

Our  rate  constant  measurements  for  the  OH  +  CH.,  H20  +  CH3  (2) 

reaction  are  plotted  along  with  selected  previous  work  in  Arrhenius  graph 

form  in  Figure  6.  Once  again  our  data  clearly  demonstrates  the  existence 

of  pronounced  Arrhenius  graph  curvature  for  Reaction  (2).  Immediately 

apparent  from  this  figure  is  the  fact  that  straight-line  extrapolations 

of  either  the  low- temperature  based  NASA  recommendation^  or  the  high- 
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temperature  based  Peeters  and  Mahnen  recommendation  produce  order  of 
magnitude  errors  in  k2(T)  in  the  opposite  temperature  regime.  Based 
principally  on  low-  and  high-temperature  experimental  measurements, 

Zellner5  derived  an  emiprical  best  fit  for  k2(T)  given  by  k2 (T )  =  2.57  x 
10~18  T2^3  exp-(1233/T)  cm3molecule-1  sec"1 .  Within  experimental  error 
both  our  data  (298-1020  K)  and  the  single  point  of  Ernst,  Wagner  and 
Zellner  (1300k)  fall  dead-on  the  curve  (solid  line  in  Figure  6)  traced 
out  by  this  expression.  Interestingly,  the  calculated  best  fit  expression 
based  on  our  data  only,  k2(T)  =  1.32  x  10-lT  T1-92  exp(-1355/T)  cm3molec-1sec_1, 
yields  a  curve  which  is  barely  distinguishable  from  that  recommended  by 
Zellner.  Indeed,  given  the  quoted  accuracies  of  these  and  previous  measure¬ 
ments  on  Reaction  (2),  a  rather  extended  array  of  parameters  could  be  sug¬ 
gested  as  producing  representative  fits  for  k2(T).  While  the  present 
experiments  have  demonstrated  both  the  appropriateness  of  the  functional 
form  k(T)  =  ATn  exp(-Eo/RT)  and  the  basic  soundness  of  previous  measurements 
in  each  separated  temperature  regime,  the  curve  fitting  calculations  indi¬ 
cate  that  extreme  caution  must  be  exercised  when  attaching  physical  significance 
to  multi-parameterized  rate  constant  expressions. 
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At  intermediate  temperatures  only  one  prior  direct  kinetic  investigation 
of  Reaction  (2)  has  been  undertaken.  Using  the  technique  of  flash  photolysis- 
resonance  absorption,  Zellner  and  Steinert  measured  k2(T)  up  to  temperatures 
of  892K.  At  lower  temperatures  these  data  agree  well  with  Zellner's 
empirical  fit,  but  at  higher  temperatures  their  measured  k2(T)  values  exceed 
the  best  fit  values  by  a  factor  increasing  with  temperature  and  already 
more  than  two  at  892K. 

One  possible  explanation  for  the  discrepancy  at  higher  temperatures  be¬ 
tween  Zellner  and  Steinert's  results  and  those  of  the  present  study  relates 
to  the  magnitude  of  the  hydroxyl  radical  reactant  concentration  demanded 
by  the  resonance  absorption  and  fluorescence  detection  modes.  Although 
both  techniques  measure  the  disappearance  of  OH  radicals  in  real  time, 
the  absorption  technique  monitors  the  difference  between  two  large  numbers 
[(It_o  -  It)/It=0]  and  thereby  requires  much  larger  OH  concentrations  to 
match  the  statistics  obtainable  from  the  fluorescence  technique.  Thus 
Zellner  and  Steinert's  initial  post-flash  OH  concentrations  were  40-100 
times  higher  and  their  [CH.*]/[0H]t_Q  ratios  were  n,  10  times  lower  than 
those  characteristic  of  the  present  study.  Such  larger  radical  concen¬ 
trations  and  lower  reactant  ratios  enhance  the  probability  of  OH  dis¬ 
appearance  due  to  processes  other  than  Reaction  (2)  and  make  more  dif¬ 
ficult  kinetic  isolation  of  this  subject  reaction.  Based  on  sensitivity 
considerations,  then,  it  appears  that  the  fluorescence  detection  tech¬ 
nique  is  best  suited  to  the  investigation  of  intermediate  temperature 
hydroxyl  radical  reaction  kinetics;  indeed,  it  is  the  authors'  opinion 
that  present  constraints  on  the  general  applicability  of  the  flash 
photolysis-resonance  fluorescence  technique  to  combustion  regime  kinetics 
of  0,  H,  and  OH  relate  more  to  maintenance  of  the  mechanical  and  optical 
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integrity  of  the  experimental  system  than  to  limitations  in  detection 
sensitivity. 

The  experiments  described  above  represent  the  first  use  of  the  flash 
photolysis-resonance  fluorescence  technique  in  the  study  of  radical- 
_ molecule  reaction  kinetics  above  500K.  The  need  for  direct  kinetic  investi¬ 
gations  of  prototype  radical-molecule  reactions  from  low  temperatures  to 
those  well  within  the  flame  region  is  indisputable.  Combustion  modeling 
requires  voluminous  rate  constant  input  from  both  direct  measurement  and 
extrapolation  [note  here  the  improvement  in  extrapolative  accuracy  achiev¬ 
able  upon  progression  of  data  availability  limits  from  400K  (1000/T  =  2.5) 
to  1000K  (1000/T  *  1.0)].  Enhancement  of  our  general  understanding  of 
bimolecular  reaction  rate  theory  is  ultimately  required  to  render  practical 
the  global  kinetic  approach  to  the  chemical  modeling  of  complex  systems. 
Unambiguous  rate  constant  studies  over  wide  temperature  ranges  (0.5  s 
1000/T  *  5.0)  will  promote  this  goal. 
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Table  I. 


Bimolecular  Rate  Constant  vs.  Temperature 
for  Reaction  {!)•  OH  +  Hg  HgO  +  H 


Temperature 

298 

499 

576 

648 

739 

838 

904 

992 


3  -1  -1 

(cmmolec  sec  ) 

1/ 

bimolecular 
(6.08  ±  0.37)  x  10*15 
(1.15  ±  0.08)  x  10~13 
(2.34  ±  0.18)  x  10'13 
(3.84  ±  0.18)  x  10~13 
(7.64  ±  0.56)  x  10‘13 
(1.30  ±  0.10)  x  10"12 
(1.86  ±  0.35)  x  10”12 
(1.99  ±  0.19)  x  10'12 


*Stated  Error  Bounds  Represent  the  2<r  Val 
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Table  II. 


Bimolecular  Rate  Constant  vs.  Temperature 
for  Reaction  (2):  OH  +  HgO  +  CH3 


3  - 1  -I 

(cm  molec  sec  ) 

k  * 

bimolecular 

(7.50  i  0.60)  x  10-15 

(4.73  4  0.45)  x  10'14 

(8.1  ±  1.1  )  x  10-14 

(1.45  4  0.12)  x  10'13 

(1.67  t  0.06)  x  10"13 

(3.14  4  0.40)  x  ID'13 

(2.75  ±  0.44)  x  10~13 

(5.78  ±  0.58)  x  10"13 

(8.4  ±  1.5  )  x  10~13 

(1.50  ±  0.15)  x  10'12 

(2.00  ±  0.20)  x  10"12 


* 


Stated  Error  Bounds  Represent  the  2 a  Values 
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FIGURE  CAPTIONS 


Fig.  1.  Schematic  drawing  of  the  experimental  apparatus.  (For  clarity 
the  flashlamp  and  resonance  lamp  are  pictured  at  180°  to  each 
other  whereas  in  reality  they  were  situated  at  90°). 

Fig.  2.  Typical  OH  temporal  profiles  observed  following  flash  photolysis 

of  H2O/CH4 /diluent  mixtures.  Experimental  conditions:  P  -  50  Torr, 
diluent  =  Ar,  flash  energy  =  120  J,  H20  pressure  =  150  mTorr. 

The  concentrations  of  methane,  temperature,  and  the  number  of 
traces  averaged  are  shown  next  to  each  line. 

Fig.  3.  Plot  of  the  first  order  rate  constant  versus  the  H2  concentra¬ 
tion  for  Reaction  (1)  at  several  temperatures. 

Fig.  4.  Plot  of  the  first  order  rate  constant  versus  the  CH4  concentra¬ 
tion  for  Reaction  (2)  at  several  temperatures. 

Fig.  5.  Arrhenius  graph  for  Reaction  (1). 

Fig.  6.  Arrhenius  graph  for  Reaction  (2). 
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The  achievement  of  an  understanding  of  the  oxidation  characteristics 
of  halogenated  boron  species  is  essential  to  the  estimation  of  the  concen¬ 
trations  of  boron  molecular  emitters  present  in  advanced  fuel  rocket  plume 
exhausts.  Knowledge  of  the  typical  concentrations  of  emitting  species  is 
a  prerequisite  to  the  interpretation  of  rocket  plume  infrared  signatures. 

Since  very  little  is  known  about  haloborane  oxidation  (and  also  about  the 
technology  required  for  its  study),  we  have  studied  relevant  molecular 
reactions  with  the  aim  of  aiding  development  of  experimental  techniques  which 
would  permit  future  investigation  of  the  more  complex  radical  reactive  systems. 

In  the  work  statement  for  this  program  we  were  contracted  to  measure 
absolute  rate  constants  for  Reactions  (1)  and  (2), 


°2  +  Bf73  Proclucts 

(1) 

C>2  +  BC1 3  -*■  Products, 

(2) 

using  the  technique  of  stop-flow  time-of-flight  mass  spectrometry.  A 
schematic  diagram  of  this  technique  is  displayed  in  Figure  1.  In  this 
system  a  thermostated  reaction  vessel  equipped  with  a  very  small  molecular 
leak  is  positioned  directly  above  the  ion  source  of  a  time-of-flight  mass 
spectrometer.  Reactants  are  admitted  to  the  vessel  and  mix  homogeneously 
and  reach  temperature  equilibrium  within  a  few  milliseconds.  Following 
mixing,  the  reaction  is  occurring  at  the  desired  conditions  of  temperature, 
pressure  and  reactor  composition.  The  molecular  leak  at  the  bottom  of  the 
reaction  cell  allows  a  very  small  fraction  of  the  cell  contents  to  diffuse 
into  the  ion  source  of  the  mass  spectrometer.  This  effluent  is  ionized,  the 
Ions  are  accelerated  to  constant  energy,  and  the  ion  current  and  time  of 
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Stop-Flow  Time-of-Flight  Mass  Spectrometer  System 


arrival  at  the  detector  are  measured.  The  entire  spectrum  can  be  cycled 
every  100  ysec.  In  essence,  this  technique  allows  the  determination  of  the 
cell  contents  in  real  time.  In  a  kinetic  experiment,  however,  only  one  ion, 
produced  upon  ionization  of  either  a  reactant  or  product  molecule,  is  monitored 

HfUh  time.  The  temperature  range  of  our  stop-flow  system  is  250K-450K. 

-23  3  -1  -1 

Second  order  recaction  rate  constants  between  2  x  10  cm  molecule  s 
and  5  x  10”17  cm3molecule"^s~^  are  amenable  to  accurate  measurement  within 
this  system.  Since  this  range  of  rate  coefficients  corresponds  to  that 
expected  for  many  molecule-molecule  reactions,  this  technique  was  proposed 
for  the  study  of  Reactions  (1)  and  (2). 

When  stop-flow  time-of-fl ight  mass  spectrometry  experiments  were  per¬ 
formed  on  Reactions  (1)  and  (2),  it  was  found  that  within  experimental  error 
addition  of  one  reactant  to  the  reaction  vessel  had  no  effect  on  the  rate  of 
disappearance  of  the  other  reactant.  We  thus  estimated  that  the  rate  con- 
stants  for  these  processes  are  -  1  x  10  cm  molecule  s”  . 

Having  found  that  the  rate  constants  for  Reactions  (1)  and  (2)  are 
too  slow  to  be  measured  by  the  above  technique,  we  devised  a  time-resolved 
infrared  spectroscopy  study  of  these  reactions.  An  optical  absorption  cell 
seasoned  with  the  haloborane  prior  to  each  experiment  was  filled  with  a  few 
torr  of  the  haloborane  and  350-500  torr  of  oxygen.  Haloborane  infrared  ab¬ 
sorption  intensities  were  then  monitored  as  functions  of  time.  Typical 
results  of  these  experiments  are  plotted  in  Figure  2.  The  haloborane  decay 
rates  are  seen  to  be  extremely  slow  and,  because  wall  loss  effects  may  have 
been  contributing  to  the  observed  decays,  the  rate  constants  extracted  from 
these  pseudo-first  order  decays  must  be  considered  as  upper  limits.  Our 
knowledge  of  these  rate  constants  thus  remains  limited  to  thi*  expressions 
ki  *  1  x  10’30  cm3molecule”^s"^  and  Y.^  -  4  x  10"^  cm3molecule"^s’^. 
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Figure  2.  Time-Resolved  Infrared  Spectroscopy 
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Having  observed  nearly  negligible  reactivities  for  the  reactions 
0^  +  BF^/BClj*  we  decided  to  probe  the  general  oxidative  behaivor  of  the 
haloboranes  in  more  detail  by  studying  their  reactions  with  more  active 
oxidizing  agents.  Investigations  of  the  following  reaction  processes  were 
Bade: 


O3  +  BF3  -k  Products 

(3) 

O3  +  BC13-  Products 

(4) 

H20  +  BC1 3  ■*  Products  . 

(5) 

Again  using  the  technique  of  time-resolved  infrared  spectroscopy.  Reaction 
(3)  was  investigated  and  found  to  be  immeasurably  slow  at  room  temperature. 

This  reaction,  then,  was  not  looked  at  further.  The  study  of  Reaction  (4), 
on  the  other  hand,  proved  to  be  very  interesting.  Time-sequenced  sets  of 
complete  infrared  spectra  of  reacting  03/BCl3/Ar  mixtures  were  taken.  The 
initial  BC1 3  and  argon  pressures  were  measured  using  a  capacitance  manometer 
and  the  initial  03  pressure  was  calculated  from  UV  absorption  measurements 
made  at  254  nm.  The  initial  reactant  partial  pressures  were  varied  by  as 
much  as  factors  of  four  in  different  runs.  The  BC1 3  and  03  decay  rates  and 
reaction  product  buildup  rate  (as  shown  in  Figure  3a,  a  reaction  product  peak 
not  observed  in  the  pure  03  or  pure  BC13  spectra  grows  with  time)  are  cal¬ 
culated  from  the  peak  intensities  in  the  sequential  infrared  spectra.  Relative 
concentration  vs.  time  after  mixing  plots  for  a  typical  experimental  run  are 
displayed  in  Figure  3b.  It  is  seen  that  the  rate  of  growth  of  the  reaction 
product  closely  images  the  rate  of  decay  of  the  BC13  reactant.  It  is  worth 
noting  here,  however,  that  the  bimolecular  rate  constants  calculated  from  such 
plots  do  show  a  slight  tendency  to  increase  with  increasing  [Oj^.q/BCIj]^ 
ratios.  These  variations  are  likely  due  to  the  overall  stoichiometry  of  the 
reaction. 
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Figure  3a.  Time-Resolved  Infrared  Spectroscopy 
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Figure  3b.  First  Order  Decay  Plot  for  the  Disappearance  of 
BC1 3  and  Appearance  of  Reaction  Product  in  the 
Reaction  Product  in  the  Reaction  of  BC1,  with 
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Three  other  experimental  techniques  were  applied  to  the  study  of 
Reaction  (4):  (1)  time-resolved  UV  absorption  spectroscopy,  (2)  stop-flow 

time-of-flight  mass  spectrometry,  and  (3)  reaction  end  product  analysis.  The 
results  of  these  studies  will  be  discussed  below  in  turn. 

A  schematic  of  the  simple  but  very  powerful  technique  of  time-resolved 
UV  absorption  spectroscopy  is  given  in  Figure  4a.  Light  from  a  mercury  lamp 
at  253.7  nm  was  passed  through  an  absorption  cell  containing  the  O^BCl^/Ar 
mixture.  The  transmitted  intensity  was  monitored  as  a  function  of  time  by 
a  filtered  photodiode.  Experiments  were  run  with  [BC1 ^3  in  varying  excesses 
over  [03],  ozone  being  the  monitored  species.  Use  of  this  technique  required 
considerable  caution.  Since  253.7  nm  radiation  is  absorbed  by  ozone  via  the 
photodissociative  reaction  0^  0(^D)  +  C^,  the  possibility  existed  that 

the  atomic  oxygen  formed  in  the  probe  process  could  contribute  to  the  chemistry 
under  study.  By  using  a  variety  of  ozone  concentrations  and  253.7  nm  radi¬ 
ation  fluxes,  however,  it  was  conclusively  shown  that  such  probe-induced 
processes  were  unimportant.  Furthermore,  by  allowing  the  low  concentration 
reactant  ozone  to  undergo  complete  reaction  and  the  measuring  essentially  no 
absorption  of  253.7  nm  radiation,  it  was  clearly  demonstrated  that  the  products 
of  the  reaction  were  not  absorbing  the  probe  radiation  and  thereby  biasing 
the  kinetic  results.  Thus  the  time  variation  of  the  transmitted  253.7  nm 
intensity  generates  an  accurate  time  history  of  the  ozone  concentration. 

The  experiments  were  run  under  pseudo-first  order  conditions  with  BCl^  in 
excess  and,  as  shown  in  Figure  4b,  the  bimolecular  rate  constant  was  obtained 

as  the  slope  of  a  plot  of  the  pseudo-first  order  decay  constant  vs.  the  BC1 3 

-21  3  -1  -1 

concentration.  A  value  of  k^.  (298K)  =  8.31  x  10  cm  molecule  s  was 
obtained. 

Reaction  (4)  was  also  studied  by  the  previously  described  technique  of 
stop-flow  time-of-flight  mass  spectrometry.  Small  concentrations  of  BCl^ 
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Figure  4a.  Time-Resolved  UV  Absorption  in  0 
Reaction. 


Figure  4b.  Plot  of  0^  Decay  Rate  vs.  [BCl^]. 
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in  varying  excesses  of  ozone  were  mixed  in  the  reaction  vessel  and  the  time 

variation  of  the  BCI^  ion  was  followed  mass  spectrometrically.  Typical  decay 

-20 

curves  are  shown  in  Figure  5a;  a  bimolecular  rate  constant  of  k^.  =  1.2  x  10 
3  -1  -1 

cm  molecule  s  was  calculated  by  the  usual  methods.  Averaging  this  result 

with  that  obtained  using  time-resolved  UV  absorption  spectroscopy  techniques, 

one  obtains  a  room  temperature  rate  constant  for  Reaction  (4)  of  3 
-20  3  -1  -1 

(1.0  ±  0.3)  x  10  cm  molecule  s  .  It  is  worth  noting  here  that  this  value 
exceeds  that  obtained  in  the  time-resolved  infrared  spectroscopy  experiments 
(see  Figure  3)  by  about  a  factor  of  six.  Though  we  do  not  feel  that  we  can 
adequately  explain  this  discrepancy,  it  is  probable  that  it  is  at  least  in 
part  related  to  wall  absorption  of  ozone  in  the  extended  duration  infrared 
spectroscopy  experiments. 

Mass  spectrometric  studies  of  Reaction  (4)  were  also  performed  by  time 
sequenced  sweeps  of  the  entire  mass  spectrum.  From  top  to  bottom.  Figure  4b 
shows  respectively  the  mass  spectrum  obtained  for  pure  ozone,  that  obtained 
for  pure  BCl^,  and  that  obtained  a  substantial  period  of  time  after  mixing 
a  small -concentration  of  BCl-j  with  a  large  excess  of  ozone.  Close  examination 
reveals  that  the  observed  product  patterns  differ  substantially  from  the 
superposition  of  the  pure  0^  and  pure  BCl^  traces.  In  Table  I  we  have  listed 
those  significant  mass  peaks  observed  in  the  mixed  O^/BCl^  system  which  were 
found  to  be  of  only  minor  importance  in  either  pure  spectrum.  Species  such 
as  B0,  B02»  BC10,  Cl2  and  HC1  have  been  identified.  Other  significant  mass 
peaks  remain  unassigned  at  this  time. 

Finally,  an  end  product  analysis  investigation  of  Reaction  (4)  was 
undertaken.  Ozone  and  BCl^  were  mixed  in  a  reaction  vessel  and  allowed  to 
react  to  completion.  The  techniques  of  infrared  and  UV  absorption  spec¬ 
troscopy  and  expansion  of  condensable  and  non-condensable  gases  into  calibrated 
volumes  were  utilized  in  the  product  analysis.  As  can  be  seen  in  Table  II, 
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Table  I. 


SIGNIFICANT  MASS  SPECTRAL  PEAKS  OBSERVED  IN  THE 
03  +  BCI3  SYSTEM  WHICH  WERE  FOUND  TO  BE  OF  ONLY 
MINOR  IMPORTANCE  IN  THE  MASS  SPECTRA  OF  PURE  03 
AND/OR  PURE  BCI 3 . 


MASS  NUMBER 

26 

27 

29 

31 

34 

36 

38 

42 

43 

51  or  52 
57 
62 

63 

64 
70 
72 
77 

no 


POSSIBLE  ASSIGNMENT 

BO 

BO 

9 

? 

9 

• 

HCI 

HCI 

B02 

bo2 

9 

9 

BCI  0 
BCI  0 
BCI  0 

Cl2 

Cl2 
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the  loss  of  BC1 ^  molecules  was  found  to  be  about  2/3  that  of  the  gain  in  CI2 
molecules:  chlorine  balance  seems  assured  through  the  formation  of  Cl2. 

Since  the  number  of  condensable  boron  oxide  molecules  formed  could  not  be 
measured,  the  boron  and  oxygen  balances  cannot  be  definitively  assigned.  One 
notices,  however,  that  the  loss  of  ozone  molecules  exceeds  the  loss  of  BC1 3 
molecules,  possibly  by  as  much  as  a  3:2  ratio.  If  such  stoichiometry  is 
assumed,  then  the  overall  Reaction  (4)  mechanism  could  be  represented  by 
reaction  expressions  of  the  type  (or  linear  combinations  thereof)  listed  at 
the  bottom  of  Table  II. 

An  abbreviated  stop-flow  time-of-flight  mass  spectrometry  study  was 
also  performed  on  Reaction  (5).  This  reaction  was  found  to  be  fast,  prob¬ 
ably  heterogeneous,  and  produced  HC1  and  condensable  boron  oxides  in  high 
yield.  Indeed,  the  addition  of  pure  BCl^  to  a  reaction  vessel  which  had 
previously  contained  H20  resulted  in  both  the  fast  decay  of  BC1 3  and  the 
obvious  formation  of  HC1.  In  experiments  in  which  BC1 3  and  H20  were  mixed 
in  the  gas  phase,  the  molecular  leak  of  the  reaction  vessel  very  quickly 
became  clogged  resulting  in  the  termination  of  gas  flow  into  the  mass 
spectrometer  analyzer. 

The  above  described  studies  of  Reactions  (1)  through  (5)  have  resulted 
in  some  interesting  information  concerning  the  general  oxidative  behavior 
of  the  haloboranes.  Equally  important,  they  have  provided  sufficient 
knowledge  to  make  possible  studies  of  more  plume- relevant  processes  such  as 
the  radical -molecule  reaction 

0(3P)  +  BF3/BC13  -*■  Products  . 
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III.  KINETIC  STUDY  OF  THE  REACTION  OF  OH  WITH 
H2  AND  D2  FROM  250  to  1050  K 


36 


INTRODUCTION 


The  high  temperature  kinetic  data  needed  for  combustion  modeling  has 
been  mainly  derived  from  shock  tube  and  flame  studies.  In  these  experiments, 
the  accuracy  of  the  obtained  rate  data  is  dependent  on  the  degree  of  under¬ 
standing  of  the  chemistry  involved  in  the  system  that  is  studied.  In  many  cases  it 
has  been  found  that  the  extrapolation  of  the  derived  high  temperature  (T  >  1000K) 
data  to  low  temperatures  (T  >  500K)  yield  rate  coefficients  which  disagree  with  the 

results  of  direct  low  temperature  measurements.  It  has  recently  been  recognized 

12  3 

that  this  discrepancy  could  be  due  to  either  non-Arrhenius  behavior  *  ’  of  the 
reaction  under  study  or  to  systematic  errors  in  the  high  temperature  data.  To 
fully  understand  the  variation  of  the  rate  coefficients  with  temperature  ir.  the 
intermediate  temperature  range,  i.e.,  298  >  T  >  1000K,  direct  experimental 
measurements  of  elementary  reaction  rate  coefficients  are  needed.  This  need 
has  recently  encouraged  adaptation  of  proven  low  temperature  kinetic  technique 
to  intermediate  and  high  temperature  studies.  The  developments  have  provided 
techniques  which  can  directly  measure  the  rate  coefficient  for  an  elementary 
reaction  over  a  wide  temperature  range,  thus  minimizing  the  chances  of 
relative  systematic  errors. 

In  parallel  to  these  experimental  developments,  great  strides  have  been 

.  4  5 

made  in  theoretical  calculations  of  elementary  reaction  rate  coefficients.  ’ 

The  transition  state  theory  has  evolved  to  a  point  where  quite  rigorous 
calculations  of  thermal  rate  coefficients  can  be  carried  out.  This  ability  can 
be  largely  attributed  to  the  availability  of  calculated  (ab  initio)  potential 
energy  surfaces  for  many  atom-diatom  reactions  and  even  four  atom  reactions 
such  as  OH  +  H2.® 

We  have  recently  adapted  the  conventional  technique  of  flash  photolysis- 
resonance  fluorescence  to  measure  OH  reaction  rate  coefficients  at  temperatures 
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7  8 

4bove  1000K.  *  The  rate  coefficient  for  the  reaction, 

OH  +  H2  +  H20  +  H  (1) 

was  measured  in  the  temperature  range  of  29S-992K  and  the  £nk^  vs.  1/T 

plot  was  found  to  be  non-linear. ^  Since  our  measurement  of  k^,  an  ab  initio 

calculation  of  the  potential  energy  surface  for  this  reaction,^  as  well  as  a 
.  .  S 

transition  state  theory  calculation  of  kj  which  utilized  this  surface,  have 
been  reported.  The  potential  energy  surface  for  the  OH  +  H2  reaction  would  be 
nearly  identical  to  that  of  OH  with  02, 

OH  +  D2  -*■  OHO  +  D  (2) 

and  therefore  comparison  of  measured  values  of  k2  with  those  calculated  would 
constitute  a  direct  test  of  the  theory.  Also,  the  question  of  how  isotopic 
substitution  affects  the  curvature  of  Arrhenius  plots  is  worth  investigating. 
Therefore,  with  the  hope  of  testing  the  theory  and  furthering  our  understand¬ 
ing  of  this  elementary  reaction,  we  have  measured  k2  over  the  temperature  range 
of  250-1050  K. 

Even  though  both  k^  and  k2  have  been  measured  by  many  investigators,  they 

have  not  been  studied  using  one  direct  technique  which  can  cover  a  wide  tempera 

ture  range.  As  mentioned  earlier,  we  have  measured  k^  in  the  temperature  range 

of  298  <  T  <  992  K.  However,  to  obtain  a  consistent  set  of  data  for  kj  and 

g 

also  since  k1  has  been  measured  only  by  one  group  at  T  <  298K,  we  extended 
the  measurement  of  k^  to  250K. 

The  technique  of  flash  photolysis-resonance  fluorescence  as  adapted  for 
high  temperature  studies  was  utilized  in  carrying  out  all  the  reported  measure¬ 
ments.  A  detailed  description  of  this  technique  is  given  in  this  paper. 
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EXPERIMENTAL  SECTION 

The  utilization  of  flash  photolysis-resonance  fluorescence  (FP-RF) 
technique  in  studies  of  OH  radical  reactions  at  temperatures  below  500K 
are  amply  described  in  the  literature. We  have  recently  adapted  this 
technique  for  studies  at  temperatures  up  to  ^  1200K.  In  this  paper  we 
will  describe  our  high  temperature  system  in  detail. 

The  primary  difficulty  in  using  the  FP-RF  technique  at  high  temperatures 
is  maintaining  good  OH  detection  sensitivity  while  using  a  reactor  which 
can  be  operated  at  these  temperatures.  We  spent  a  considerable  amount  of 
time  experimenting  with  various  cell  designs,  cell  materials,  focusing  optics, 
and  light  baffle  arrangements  to  maximize  OH  detection  sensitivity  and  still 
be  able  to  heat  the  reactor  to  high  temperatures.  The  final  apparatus  con¬ 
figuration  that  was  adopted  is  shown  in  Figure  1.  The  principal  system 
components  are  (1)  a  thermostated  reaction  cell,  (2)  a  spark  discharge  flash- 
lamp  perpendicular  to  one  face  of  the  cell,  (3)  a  CW  OH  resonance  lamp 
perpendicular  to  the  photolysis  beam,  (4)  a  bandpass  filter/photomultiplier 
combination  for  monitoring  OH  resonance  fluorescence  emitted  perpendicular 
to  both  the  photolysis  and  resonance  radiation  beams,  and  (5)  a  signal 
averager  and  fast  photon  counting  electronics. 

The  reactor  was  built  out  of  quartz.  Five  highly  polished  Suprasil  windows 
were  fused  to  the  cell  using  quartz  tape.  This  cell  was  coated  with  graphite 
on  the  outer  surface.  Once  baked,  the  graphite  bonded  well  with  the  quartz 
surface.  This  black  graphite  coating  served  two  purposes:  (1)  it  cut 
down  the  amount  of  scattered  radiation,  and  (2)  it  formed  a  good  thermally 
conductive  surface  which  helped  avoid  'hot  spots'.  This  cell  was  heated 
using  tantalum  wire  windings  which  were  contained  inside  ceramic  tubes.  The 
heaters  were  powered  by  *  constant  voltage  power  supply.  The  temperature 
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could  be  controlled  by  changing  the  applied  voltage  using  variacs.  The 
cell  and  the  insulated  heaters  were  covered  with  at  least  5  layers  of  thin  stain 
less  steel  sheets  which  acted  as  radiation  shields,  thereby  minimizing  the  heat 
loss  from  the  cell.  A  chromel-alumel  thermocouple  was  attached  to  the  outer 
.surface  of  the  reactor  primarily  to  aid  in  knowing  the  temperature  of  the  cell 
during  heating  and  cooling  periods.  The  temperature  of  the  gas  inside  the 
reactor  was  directly  measured  using  a  thermocouple  which  was  inserted 
into  the  reactor.  In  preliminary  characterization  experiments  the 
temperature  differential  between  the  upper  and  lower  portions  of  the 
reaction  zone  (maximum  length  ^  2  cm)  was  found  to  be  <  5K  at  1000K  and 
unmeasurable  at  500K.  The  bottom  window  (fluorescence  port  located  ^  5  cm 
below  the  active  volume)  was  measured  to  be  ^  15K  below  that  of  the  active 
volume  at  1000K  when  the  gases  were  flowing  through  the  cell.  Therefore,  it 
can  be  quite  confidently  ascertained  that  at  1000K  the  temperature  in  the 
active  zone  was  constant  to  within  5K  and  known  to  better  than  1  OK. 

Quartz  (unlike  pyrex)  transmits  310  nm  radiation  and  thus  increases  the 
background  scattered  resonance  lamp  radiation  that  reaches  the  detector.  To 
overcome  this  problem,  two  modifications  were  carried  out.  (a)  The  microwave 
discharge  lamp,  the  source  of  resonance  radiation,  was  moved  away  (^  80  cm) 
from  the  cell.  A  5  cm  diameter  quartz  convex  lens  was  placed  20  cm  from 
the  discharge  to  collect  the  radiation  and  render  it  nearly  parallel.  A 
second  5  cm  diameter  quartz  lens  located  n,  15  cm  from  the  center  of  the 
reactor,  focused  this  parallel  beam  on  the  far  port  of  the  reactor.  This 
configuration  enabled  almost  all  reasonance  radiation  to  pass  only  through 
the  cell  windows,  (b)  The  photon  counting  electronics  were  improved  to  count 
linearly  up  to  3  MHz.  With  this  combination,  the  measured  background  count 
rate  ranged  from  100-200  KHz.  (It  should  be  pointed  out  that  generally 
emphasis  is  placed  on  reducing  the  background  rate  to  a  very  low  value,  often 


at  the  expense  of  decreasing  the  amount  of  resonance  radiation  that  enters 
the  reactor.)  The  output  of  the  N2  flashlamp  was  collimated  by  a  set  of 
baffles  and  then  rendered  parallel  by  a  pith  polished  Suprasil  lens  which 
ensured  that  all  flash  radiation  exited  only  through  the  windows.  In  addi¬ 
tion,  a  set  of  lenses  was  employed  to  collect  the  OH  fluorescence  from  the 
center  of  the  cell  and  then  focus  it  onto  the  photocathode  of  a  RCA  8850 
photomultiplier  tube.  Use  of  this  lens  set  resulted  in  increases  in 
both  fluorescence  signal  and  background  light  with  a  large  improvement  in 
signal-to-noise  ratio. 

To  determine  the  sensitivity  of  the  apparatus,  we  carried  out  back-to- 
back  experiments  at  298K  where  the  signal/noise  per  flash  was  determined  under 
typical  experimental  conditions  (flash  energy  =  70J,  50  mTorr  H20,  100  Torr  Ar) 
and  compared  with  the  signal/noise  per  flash  obtained  under  conditions  where 
a  known  OH  concentration  could  be  produced.  (Here,  the  signal/noise  is  defined 
as  the  ratio  of  signal  due  to  OH  to  the  square  root  of  the  background  counts.) 
The  procedure  used  to  generate  a  known  concentration  of  OH  involved  replacing 
the  flashlamp  with  a  frequency  quadrupled  pulsed  266  nm  laser  beam  (from  a  Nd: 
YAG  laser)  and  adding  ^  5  mTorr  0^  to  the  reaction  mixture.  Upon  photolysis 
the  following  reactions  take  place: 


°3  26£) 


hv 


nm 


*  0( 1D)  +  02(]A  ) 


I  0(3P)  +  o2(3e) 


O^D)  +  Ar  -  0(3P)  +  Ar 
O^D)  +  H20  -*•  2  OH  . 


The  laser  beam  was  expanded  so  its  cross-sectional  area  as  it  traversed  the 
reaction  cell  was  approximately  equal  to  that  of  the  flashlamp  radiation. 

The  laser  power  was  measured  at  the  rear  of  the  cell  with  an  absolute  energy 
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monitor  and  corrected  for  the  loss  due  to  reflection  by  the  back  window 

1 2 

(%  10',).  Since  the  absorption  cross  section  for  0^  at  266  nm,  the  quanto" 
yield  for  O^D)  from  03  at  266  nir.  0.9),]3  the  rate  coefficients  for  0(1D) 
quenching  by  Ar^  and  reaction  with  ^0^  are  known,  we  could  determine  the 
concentration  of  OH  in  the  reactor. 

The  results  of  actinometry  indicate  that,  for  the  reaction  mixture 
employed,  %  1  x  109  OH  per  cm3  could  be  detected  (signal/noise  *  1)  with  an 

_3 

integration  time  of  1  x  10  s.  Extrapolation  of  this  result  to  other  experi¬ 
mental  conditions  shows  that  all  experiments  were  carried  out  with  OH  con- 

10  11  *> 

centrations  between  1  x  10  and  1  x  10  per  cm''.  It  was  also  found  that  as 
long  as  the  reactor  was  kept  clean,  the  measured  signal  for  a  given  OH  con¬ 
centration  did  not  deteriorate  significantly  with  increases  in  temperature. 

All  experiments  were  carried  out  under  pseuto-first  order  kinetic  condi¬ 
tions  with  [H£]  and  [D2]  in  large  excess  over  [OH].  (^J/COH]  and  [D2]/[0H] 

ratios  were  always  greater  than  1000.  [H2]  ranged  f rom  ^  2  x  1013  to  2  x  1017  cm' 

13  18  -3 

and  [D2]  ranged  from  ^  2  x  10  to  1  x  10  cm  .)  OH  radicals  were  produced  by 
flash  photolysis  of  H20  at  wavelengths  between  the  Suprasil  cutoff  at  165  nm  and 
the  onset  of  continuum  absorption  at  185  nm  (flash  duration  *  50  psec). 

Following  the  flash,  weakly  focused  OH  resonance  lamp  radiation  continuously 
excited  a  small  fraction  of  the  OH  in  the  reactor  to  the  electronically 
excited  A  Z  state*,  ths  resultant  (0,0)  band  A  -*■  X  fluorescence  emanating 
in  the  direction  perpendicular  to  both  the  resonance  excitation  and  photolysis 
beams  was  collected  by  a  lens  and  focused  onto  a  photomultiplier  fitted  with 
a  bandpass  filter  (309.5  nm  peak  transmission,  10  nm  FWHM).  Signals  were 
obtained  by  photon  counting  and  then  fed  into  a  signal  averager  operated  in 
the  multichannel  scaling  mode.  For  each  decay  rate  measured,  sufficient 
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flashes  (50-1000)  were  averaged  to  construct  a  well-defined  temporal  profile 
over  at  least  a  factor  of  twenty  variation  in  [OH], 

All  experiments  were  carried  out  under  "slow  flow"  conditions.  The  flow 
rate  through  the  reactor  was  such  that  each  photolysis  flash  encountered  a 
fresh  reaction  mixture  (photolysis  repetition  rate  n,  0.3  Hz).  At  temperatures 
greater  than  270K,  82^)  was  flowed  from  a  12S.  source  bulb  containing  ( D 2)/ 
diluent  gas  at  800  Torr  through  distilled  water  at  room  temperature.  The 
H2(D2)  diluent  gas  mixture,  H20/diluent  gas  mixture,  and  additional  diluent 
gas  were  mixed  prior  to  entering  the  reaction  cell.  Concentration  of  each 
component  in  the  reaction  mixture  was  determined  from  measurements  of  the 
appropriate  mass  flow  rates  (measured  using  calibrated  mass  flow  meters)  and 
the  total  pressure.  At  temperatures  below  270K,  the  amounts  of  and 
needed  to  carry  out  the  experiments  were  so  large  that  pure  H2  or  02  was 
needed  in  the  source  bulb.  Therefore,  we  prepared  large  amounts  of  H2(D2)/ 
diluent  mixtures  in  an  ^  4(K  volume  and  flowed  it  directly  into  the  reactor 
after  mixing  it  with  a  small  amount  of  H20/diluent  mixture.  This  procedure 
avoided  any  possible  flow  meter  calibration  errors. 

The  gases  used  in  this  study  had  the  following  stated  purities: 

Ar  >  99.9995*,  H2  >  99.9999*,  02  >  99.96*.  Also,  the  mixture  had  less  than 
1  ppm  of  organics  (mainly  as  CH^)  and  an  isotopic  purity  of  99.5*.  and 
02  were  passed  through  liquid  nitrogen  traps  filled  with  glass  beads  before 
being  used  to  remove  any  condensible  organics.  Ar  was  used  as  supplied. 

RESULTS 

All  experiments  were  carried  out  under  pseudo-first  order  conditions  with 
[H2]  or  [D2]  »  [OH].  The  temporal  profile  of  [OH],  then,  is  given  by  the  equation, 

£n{[0H]o/[0H]t)  =  (kx[x]  +  kd)t  =  k't 
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where  X  is  either  or  D2,  k^  is  the  first  order  rate  constant  for  the  loss 

of  OH  in  the  absence  of  X,  and  k'  is  the  measured  pseudo-first  order  rate 
constant.  [OH]  decays  were  observed  to  be  exponential  for  at  least  three 
1/e  times  as  typified  by  data  shown  in  Figure  2;  the  exponential  [OH]  decay 
validates  the  existence  of  first  order  conditions.  At  each  temperature,  the 
[OH]  decay  rates  were  measured  at  several  (minimum  4,  on  the  average  6)  fixed 
reactant  concentrations.  The  bimolecular  reaction  rate  constant  was  obtained 
from  the  slope  of  the  k'  vs.  reactant  concentration  plots.  Figure  3  shows 
one  such  plot  for  k^  at  932K.  The  measured  values  of  kj  and  k 2  at  various 
temperatures  are  listed  in  Tables  I  and  II,  respectively.  The  quoted  errors 
are  2a  and  represent  the  precision  of  the  measured  rate  coefficient.  It  is 
estimated  that  the  accuracy  at  298K  is  ^  122  decreasing  to  ^  20 2  as  the 
temperature  is  either  raised  to  1000K  or  lowered  to  250K. 

It  has  been  previously  shown  that  under  our  experimental  conditions 
radical-radical  reactions  such  as  H  +  OH  or  OH  +  OH  are  insignificant.11  To 
ensure  that  our  present  results  were  devoid  of  such  secondary  reactions 
the  [0H]q,  the  flash  energy,  the  [H20],  and  the  system  pressure,  were 
varied  at  most  temperatures.  The  results  indicated  the  absence  of  any  secondary 
reactions  contributing  to  the  measured  value  of  k-j  and  k^. 

The  values  of  k-|  and  k2  are  quite  small  at  low  temperatures.  Therefore, 

-4 

very  small  concentrations  of  reactive  impurities,  in  the  order  of  10  %,  can 
contribute  significantly  to  the  measured  values  if  they  are  present.  It  is 
always  very  hard  to  completely  rule  out  this  possibility  of  impurity  reactions. 
Our  samples  of  H2  and  D2  are  very  pure.  Both  H2  and  D2  have  less  than  1  ppm 
of  CH^,  the  only  organic  that  was  present  in  measurable  quantities.  Great 
care  was  taken  to  ensure  that  no  impurities  were  added  in  the  course  of 
transferring  H2  and  D2  from  the  tanks  to  the  reactor.  Therefore  we  believe 
that  our  values  of  and  k2  are  accurate  to  at  least  202  at  250K. 
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Figure  4  shows  plots  of  £nk^  and  £nk^  vs.  1/T.  Results  of  our  previous 
measurements  of  k-j  in  the  temperature  range  of  298  to  1000  K  are  also 
included.  It  is  clear  that  Arrhenius  plots  are  not  linear  for  both  k^ 
and  k£;  tnk^  vs.  1/T  being  more  nearly  linear  than  Jtnk-j  vs.  1/T.  We  pre¬ 
viously  fit  k-j  values  between  298K-1000K  to  the  function  k(T)  =  4.12  x  10"^ 
T2'44  exp(-1281/T)  while  Zellner2  obtained  k(T)  =  1.66  x  10'16  T1,6  exp(-1660/T) 
based  on  all  previous  measurements.  Our  present  data  is  best  described  by 
the  former  of  the  two  expressions.  A  similar  three-parameter  expression 
for  Reaction  (2)  based  solely  on  our  data  is,  k^(T)  =  4.37  x  10  ^  T^‘^ 
exp(-2332/T). 
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DISCUSSION 

Reaction  (1)  has  been  studied  by  numerous  investioators"'^',"'^’''^,''°  at 

room  temperature.  There  are  four  studies  of  as  a  function  of  tempera- 

9  20  °1  22  21 

ture.  ,c  ’  With  the  exception  of  Westenberg  and  de'riaas,  all  previous 

studies  had  an  upper  temperature  limit  of  500K.  Table  ijl  lists  both  the 

29SK  value  and,  where  appropriate,  the  Arrhenius  parameters.  There  is 

reasonable  agreement  between  various  (293K)  values.  As  pointed  out  earlier  the 

£nk-j  vs.  1/T  (Arrhenius)  plot  is  not  linear.  For  the  sake  of  comparison  we  have 

fitted  our  low  temperature  (T  <  400K)  data  to  an  Arrhenius  expression  which 


yields; 


k(T)  =  (4.9  ±  o.5)  x  10-12  exp{ (1 . 99  ±  0.34)  x  103/T)  cm3molecule"  V1 


where  the  errors  are  2o  and  =  A  £nft.  These  values  are  listed  in  Table  111. 

Our  value  of  A  and  E/R  are  in  good  agreement  with  other  measurements  except 

g 

with  Smith  and  Zellner's  values.  Smith  and  Zellner  were  the  only  investigators 
who  have  measured  k^  at  T  <  298K.  They  encountered  some  difficulty  at  tempera¬ 
tures  below  270K  where  they  had  secondary  reactions  such  as 


OH  +  OH  +  M  -*•  H202  +  M  (6) 

and 

OH  +  H  +  M  -  H20  +  M  .  (7) 

Contributing  to  the  measured  OH  temporal  profile.  They  corrected  their 

measured  value  to  take  into  consideration  the  contributions  of  these 

secondary  reactions.  The  value  of  2kg  +  k^  they  used  to  correct  their 

data  was  n,  7  x  10'3^  cm®molecule”2s”^  which  in  the  light  of  recent  measure- 
23 

ments  seems  rather  high.  Therefore  it  can  be  postulated  that  they  over 

corrected  their  low  temperature  value  which  could  account  for  the  larger 

20 

E/R  values.  Figure  4  shows  the  Arrhenius  plots  reported  by  Greiner, 
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9  22 

Smi th  and  Zellner,  and  Atkinson  et  al.  It  is  clear  from  Figure  4  tnat  our 
individual  data  points  at  T  *  29£K  are  in  excellent  agreement  with  the  line  of 
Smith  and  Zellner.  It  is,  however,  our  data  T  >  29SK  which  are  lower  than  that 
of  Snith  and  Zellner's  line.  Therefore  it  is  doubtful  that  the  low  temperature 
corrections  of  Snith  and  Zellner  could  account  for  the  discrepancy.  Also,  at 
T  29SK,  the  results  of  Kestenberg  and  deKaas,  Atkinson  et  al . ,  and  Greiner 
are  all  in  good  agreement  with  our  values  of  k..  It  is  particularly  interesting 
to  note  the  agreement  over  the  entire  temperature  range  of  our  results  with 
that  of  Westenberg  and  deKaas  who  were  the  first  to  notice  the  non-linear  varia¬ 
tion  of  £nk^  with  1/T. 

18  24  25 

Reaction  (2)  has  been  measured  by  three  separate  groups  all  at  298K.  *  * 

0 

The  temperature  dependence  of  k2  has  been  studied  only  by  Smith  and  Zellner. 

Table  IV  lists  all  available  data  on  Reaction  (2).  Again  for  the  sake  of  com¬ 
parison,  we  have  fit  our  low  temperature  data  to  an  Arrhenius  form  and  obtained 
A  and  E/R  which  are  shown  in  Table  IV.  There  is  reasonable  agreement  between  all 
measured  values  of  k2  (298K);  our  value  being  the  lowest.  Even  the  Arrhenius 
parameters  derived  in  this  study  based  on  low  temperature  (i.e.  T  <  47DK)  data, 

k2(T)  “  (1.21  ±  0.52)  x  10"12  exp{  (-2.67  ±  .15)  x  103/T)  cm3molecu  e'V1  agree 

g 

well  with  that  of  Smith  and  Zellner. 

One  of  the  main  aims  of  this  work  was  to  obtain  a  set  of  data  using  a 
single  technique  over  a  wide  temperature  range  for  a  set  of  related  reactions 
such  that  a  coherent  data  base  will  be  available  to  compare  present  day  theory 
with  experimental  results.  Based  on  our  data  two  observations  can  be  made. 

(a)  The  £nk  vs.  1/T  plots  are  not  linear  for  both  Reactions  (1)  and  (2).  (b)  The 

kinetic  isotopic  effect,  i.e.,  kj/k2  increases  with  decreases  in  temperature. 
Comparison  of  experimental  results  with  theory  can  be  made  at  different  levels 
of  quantification.  The  general  non-linear  Arrhenius  behavior  for  reactions 
such  as  Reactions  (1)  and  (2)  can  be  qualitatively  explained,  and  is  even  to  be 
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expected,  based  or,  chemical  kinetic  theories  of  elementary  reactions.  Zellner 
3 

and  Gardiner  have  discussed  the  probable  causes  for  this  behavior  which  car.  be 

understood  tc  different  levels  of  precision  based  on  the  transition  state  theory 

and/or  the  collision  theory. 

Before  comparing  the  theoretical  calculations  with  experimental  results, 
the  mechanism  for  the  reaction  of  OH  with  hydrogen  needs  to  be  explicitly  stated. 

This  reaction  seems  to  proceed  through  the  bond  scission  in  hydrogen. 


HO  +  H.  -  H.  -  H  OH.  +  H. 
a  b  b  a  b  b 


rather  than  the  0-H  bond  cleavage 


+  H.  -  HK  -  H.  -  0  -  H.  +  Ha  . 
a  b  b  b  b  a 


This  assumption  is  made  in  all  calculations  and  is  seldom  stated  explicitly. 

There  is  evidence  to  justify  this  mechanism  based  on  the  work  of  Westenberg 
21 

and  deHaas  who  observed  the  almost  exclusive  production  of  H  atoms  in  the 
reaction  of  OD  with  (rather  than  D  atoms  which  would  be  produced  if  the 

25 

bond  in  hydroxyl  radical  was  cleaved)  and  the  work  of  Paraskevopoulos  and  Nip, 
who  noted  the  absence  of  kinetic  isotopic  effect  between  the  reaction  of  OH 
with  H^  and  OD  with 

The  transition  state  theory  can  be  used  at  many  levels  of  rigor  to  calcu¬ 
late  the  value  of  k^(T)  and  k2 (T)  (and  hence  (T)/k2(T)).  Assuming  a  transition 
state  structure  one  can  predict  the  temperature  variation  of  the  bimolecular 
rate  coefficient  k,  especially  when  the  transition  state  properties  are  made  to 
fit  the  value  of  k  measured  at  one  temperature.  Recently,  however,  ab  initio 
potential  energy  surfaces  have  been  calculated  for  simple  atom-diatom  reactions 
and  even  four  atom  reactions  such  as  Reactions  (1)  and  (2).®  Using  this 

potential  energy  surface,  conventional  transition  state  theory  calculations 

5  26 

have  been  carried  out  by  Schatz  and  Walch  and  Isaacson  and  Truhlar  for 

71  26 

Reaction  (1)  and  by  Schatz  and  Wagner  and  Isaacson  and  Truhlar  for  Reaction  (2). 

26 

Isaacson  and  Truhlar  have  also  calculated  k-j(T)  and  k2(T)  based  on  canonical 


variation  theory.  Figures  (5)  and  (6)  show  the  calculated  values  of  k-j(T) 
and  k^T),  respecti vely,  based  on  conventional  transition  state  theory  which 
includes  and  excludes  tunneling  as  well  as  the  canonical  variation  theory 
which  excludes  tunneling.  Two  points  are  clear:  (a)  calculated  k^(T)  values 
are  much  lower  than  experimental  values  when  tunneling  is  neglected,  the  dis¬ 
crepancy  between  calculated  and  experimental  values  increasing  with  decreasing 
temperature,  and  (b)  the  k2 (T)  values  which  are  calculated  on  the  same  potential 

energy  surface  as  that  for  Reaction  (1)  and  which  includes  tunneling  correction, 

27 

predicts  k2 (T )  values  much  lower  than  measured  values.  The  fact  that  the 

calculated  rate  coefficients  are  very  sensitive  to  the  exact  choice  of  input 

5 

parameters  is  seen  when  we  compare  the  results  of  Schatz  and  Kalch  with 

26 

those  of  Isaacson  and  Truhlar.  Using  the  same  potential  energy  surface, 
they  obtain  different  results  which  are  directly  attributed  to  differences 
in  partition  functions,  and  the  incorporation  of  anharmonicity,  spin  orbit 
coupling,  etc.  The  canonical  variation  theory  calculations  which  exclude 
tunneling  yield  the  lowest  calculated  values  for  both  k^(T)  and  k^T).  In 
fact,  the  difference  between  the  calculated  and  measured  values  of  k^(T) 
at  250K  is  ^  a  factor  of  20.  However,  it  is  expected  that  inclusion  of 
tunneling  would  greatly  enhance  the  agreement  between  calculations  and 
experimental  results.  It  should  be  noted  that  all  these  calculations  suggest 
that  tunneling  is  a  major  contributor  to  the  overall  rate  coefficient  for  the 
reaction  and  the  tunneling  contribution  increases  quite  dramatically  with 
decreases  in  temperature. 

Another  experimental  parameter  which  can  be  compared  with  theoretical 
calculations  is  the  kinetic  isotope  effect,  i.e.,  k1/k2.  Table  V  lists  these 
values  for  all  the  cases  considered  above  as  well  as  the  calculations  of  Smith 

Q 

and  Zellner.  Smith  and  Zellner  used  two  different  semiempirical  surfaces. 
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i . e - ,  LEPS  and  6EB0,  and  they  neglected  tunneling.  It  should  also  be  noted 
that  they  assured  the  energy  barrier  height  to  be  the  sare  as  the  experi¬ 
mentally  measured  activation  energy.  It  is  clearly  seen  that  all  calculations 
show  the  observed  trend  in  k^/k^,  i.e.,  decrease  with  increase  in  temperature. 
Moreover,  it  is  interesting  to  see  that  the  best  agreement  between  measured 
and  calculated  values  of  k-j /k^  is  obtained  when  tunneling  is  neglected.  Even 
in  the  case  of  the  semiempirical  calculation  of  Smith  and  Zellner,  there  is 
good  agreement  between  measured  and  calculated  kinetic  isotope  effects.  In 
fact,  this  surprising  agreement  led  to  the  speculation  that  tunneling  is 
probably  not  important  in  Reactions  (1)  and  (2).  However,  in  spite  of  the 
good  agreement  in  isotope  effect  ratios,  the  absolute  values  of  k-j  and  k2 
are  grossly  underestimated  when  tunneling  is  neglected  as  shown  by  calculations 
of  Schatz  et  al.^’^  and  Truhlar  et  al 

At  this  point  one  is  forced  to  ask  the  question,  even  if  an  accurate  ab. 
ini tio  potential  energy  surface  is  available,  can  the  present  theories  accurately 
predict  thermal  rate  coefficients,  and,  can  the  theories  shed  more  light  on  the 
exact  nature  of  the  reactive  process.  Based  on  the  calculations  that  are 
presented  it  is  seen  that  even  though  great  strides  have  been  made  in  this 
direction,  the  theory  still  needs  to  be  refined  to  predict  rate  coefficients 
with  the  accuracy  of  1 0-50%.  Of  course,  it  should  always  be  kept  in  mind  that 
ultimately  the  agreement  between  theory  and  experiment  is  a  reflection  on  the 
accuracy  of  the  potential  energy  surface  itself. 

The  explanation  of  the  non-Arrhenius  behavior  of  k^  and  k2  based  on  the 

2  3 

collision  theory  has  been  discussed  by  Zellner  and  Gardiner.  As  pointed  out 
by  Zellner,  if  the  rate  coefficient  for  the  reaction  of  0H(v=l)  with  H2  is  'v 
150  times  that  for  H2(v*0),  then  depending  on  the  efficiency  with  which  the 
Increased  vibrational  energy  can  help  hasten  the  reaction,  the  observed 
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thermal  rate  coefficient  can  exhibit  different  degrees  of  non-Arrhenius 

28 

curvature.  Zellner  et  al .  have  observed  that  H2(v=l)  is  renoved  by  OH 

-12  3  -1  -1 

with  a  rate  coefficient  of  n,  1  x  10  cm  molecule  s  .  If  this  removal  is 
attributed  totally  to  reaction  rather  than  vibrational  relaxation,  a  pos¬ 
sibility  that  is  not  experimentally  proven,  then  H2(v=l)  can  contribute 
significantly  to  k-j (T }  even  at  temperatures  as  low  as  800K.  Of  course, 
this  mechanism  does  not  preclude  contribution  due  to  tunneling  at  low 
temperatures.  In  the  absence  of  any  data  on  D2 ( v=l )  it  is  hard  to  speculate 
if  such  a  process  can  contribute  to  k^.  It  should  be  noted  that  the  pop¬ 
ulation  of  in  v=l  would  be  much  higher  than  that  of  H2  and  if  D2(v=l) 
reacts  as  fast  as  H2(v=l)  is  postulated  to  do,  the  contribution  of  D2(v=l) 
can  be  more  than  that  of  D2(v=0)  at  temperatures  as  low  as  900K,  and  hence 
would  be  easier  to  measure. 

Finally,  it  is  worth  inquiring  what  experimental  work  can  be  carried  out 
to  further  check  and  aid  in  the  development  of  the  theories  of  elementary 
reactions.  Two  experiments  are  suggested,  (a)  extension  of  measurement  of 
k-j  and  k2  to  temperatures  much  lower  than  250K  where  tunneling  should  dominate 
and  temperatures  higher  than  1000K  where  energetics  are  less  important  and 
(b)  measure  rate  coefficients  for  the  reaction  of  OH  with  HD.  Since  the 
same  potential  energy  surface  that  is  used  for  H2  and  D2  will  be  applicable, 
it  will  be  interesting  to  see  if  the  theory  can  predict  the  measured  thermal 
rate  coefficients  for  the  reaction  of  OH  with  HO.  Both  these  experiments 
will  be  underway  shortly  in  our  laboratory. 
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Table  I.  Rate  Coefficient  vs.  Temperature 
for  Reaction  1:  OH+f^-^O+H 


Temperature 

(K) 

3  -1  -1 

,  cm  molecule  s 

250 

(1.69  ±  0.14)  x  10"15 

265 

(2.67  ±  0.37)  x  10"15 

282 

(4.14  ±  0.30)  x  10'15 

295 

(5.64  ±  0.60)  x  10"15 

333 

(1.28  ±  0.05)  x  10‘14 

356 

(1.83  ±  0.08)  x  10"14 

366 

(2.02  ±  0.21)  x  10‘14 

396 

(3.32  ±  0.11)  x  10'14 

960 

(2.40  ±  0.25)  x  10-12 

1050 

(3.20  ±  0.44)  x  10‘12 
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Table  II.  Rate  Coefficient  vs.  Temperature 
for  Reaction  2:  OH+D^-’-OHO+H 


Temperature 

(10 

3  -1  -1 

k^,  cm  molecule  s 

250 

(3.16  ±  0.25)  x  lO*16 

268 

(5.00  ±  0.57)  x  10'16 

298 

(1.83  ±  0.12)  x  10'15 

336 

(3.98  ±  0.17)  x  10'15 

356 

(6.03  ±  0.36)  x  10'15 

366 

(7.07  ±  0.71)  x  10-15 

380 

(1.07  i  0.05)  x  ID"14 

398 

(1.18  ±  0.18)  x  10-14 

438 

(3.07  ±  0.23)  x  10'14 

449 

(3.35  ±  0.14)  x  10'14 

471 

(4.81  ±  0.37)  x  10-14 

525 

(8.86  ±  0.85)  x  10"14 

572 

(1.74  ±  0.06)  x  10"13 

605 

(1.61  ±  0.10)  x  10"13 

617 

(2.31  ±  0.04)  x  10'13 

664 

(3.47  ±  0.20)  x  10'13 

752 

(4.82  ±  0.21)  x  10'13 

788 

(5.93  ±  0.34)  x  10"13 

820 

(7.8  ±  1.1)  x  10'13 

932 

(1.11  ±  0.03)  x  10"12 

960 

(1.28  ±  0.11)  x  10-12 

1050 

(1.61  ±  0.45)  x  10"12 
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Table  III.  Summary  of  Rate  Constant  Data  for  Reaction  1;  0H+H.-*-H90+H 
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Observed  curved  tnk  vs.  1/T  plots,  and  hence  did  not  report  E/R  values. 
This  expression  has  been  derived  from  only  the  low  temperature  data. 


Table  IV.  Summary  of  Rate  Constant  Data  for  Reaction 


Table  V.  Comparison  of  Experimentally  Measured  Kinetic  Isotope  Effect 
(K-i  (T)/k?(T) )  with  Calculated  Values. 


Figure  Captions 


1.  Schematic  diagram  of  the  high  temperature  flash  photolysis-resonance 
fluorescence  apparatus.  AD.  amplifier  discriminator;  D,  diluent  gas; 
DVM,  digital  voltmeter;  FI,  309.5  nm  bandpass  filter;  FL  flashlamp; 

FT,  flow  transduder;  HC,  high  voltage  capacitor;  MC,  mixing  chamber; 
MCA,  multichannel  analyzer;  MG,  microwave  generator;  NV,  needle  valve; 
PG,  pressure  gauge;  PM1 ,  photomultiplier  RCA  8850;  PS,  high  voltage 
power  supply;  R/D,  reactant/diluent  gas;  RL  resonance  lamp;  TC, 
thermocouple;  TTY,  teletype;  VM,  vacuum  housing. 

2.  Typical  [OH]  temporal  profile  following  flash  photolysis  of  H^O/D^/Ar 

mixture.  Experimental  conditions:  T  =  932K;  P  =  100  Torr  (Ar); 

Flash  Energy  =  80J;  [H^O]  =  70  mTorr;  Concentration  of  in  molecules 

cm'3,  (a)  1.08  x  1014,  (b)  2.09  x  1014,  (c)  4.88  x  1014. 

3.  Plot  of  the  pseudo-first  order  rate  constant  vs.  [D2]  at  932K.  Solid 
line  was  obtained  by  linear  least  squares  analysis  of  k'  vs.  [D?]  date. 
Flash  energy,©  -  80J,  □-  200J. 

4.  Plots  of  tnk^  and  £, nk^  vs.  1000/T.  For  the  sake  of  clarity,  individual 

error  bounds  are  not  shown,  o  -  k^  values  from  Ref  (7);  □  -  k^  values 

from  Ref  (21  \  Arrhenius  parameters  reported  by  previous  investigators 
are  as  follows:  — ,  Smith  and  Zellner  Ref  (9);  —  Atkinson  et  al . 

Ref  (22 X  —  Greiner  Ref  (20). 

5.  Comparison  of  experimental  values  of  k^(T)  with  calculations  carried 

out  using  the  ab  initio  surface  of  Dunning  et  al .  Ref  (6). 

—  TST  calculation  including  tunneling  correction  (Schatz  and  Walch 
Ref.  5);  ---  TST  calculation  excluding  tunneling  (Schatz  and  Walch 
Ref.  5);  — x—  TST  calculation  including  tunneling  correction 

(Isaacson  and  Truhlar  Ref.  26);  -  TST  calculation  excluding 

tunneling  (Isaacson  and  Truhlar  Ref.  26);  — c—  Canonical  variation 
theory  calculation  excluding  tunneling  (Isaacson  and  Truhlar  Ref.  26). 

6.  Comparison  of  experimental  values  of  k^(T)  with  calculations  carried 

out  using  the  ab  initio  surface  of  Dunning  et  al .  Ref  (6).  —  TST 

calculation  including  tunneling  correction  (Schatz  and  Wagner  Ref.  27); 
—x—  TST  calculation  including  tunneling  correction  (Isaacson  and 

Truhlar  Ref.  26);  -  TST  calculation  excluding  tunneling  (Isaacson 

and  Truhlar  Ref.  26);  — c—  Canonical  variation  theory  calculation 
excluding  tunneling  (Isaacson  and  Truhlar  Ref.  26). 
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KINETICS  OF  THE  OH  REACTIONS  WITH 
CO,  HC1 ,  C2Hg  AND  C3Hg 


The  reaction  of  hydroxyl  radicals  with  CO, 

OH  +  CO  -  C02  +  H  (1) 

is  of  extreme  importance  in  combustion  and  atmospheric  processes.  In  com¬ 
bustion  systems.  Reaction  (1)  is  the  principal  mode  of  conversion  of  CO  to 
C02.  Since  this  reaction  is  very  exothermic,  its  kinetics  considerably 
impact  the  time  profile  for  combustive  energy  release,  and  thereby,  factors 
such  as  ignition,  acceleration  and  propogation  of  flames.  Also,  in  the 
atmosphere.  Reaction  (1)  is  the  principal  CO-*-  C02  conversion  pathway;  its 
kinetics  thus  impact  numerous  important  atmospheric  chemical  cycles  (e.g., 

NO  -  N02). 

The  reaction 

OH  +  HC1  -  H20  +  Cl  (2) 

is  one  of  the  major  pathways  by  which  HC1  is  destroyed  in  plumes  and  in  the 


stratosphere. 

The  reactions 

OH  +  C2Hfi  -  C2H5  +  H20 

(3) 

and 

OH  +  C3H8  C3H7  +  H20 

(4) 

are  important  in  hydrocarbon  oxidation. 

As  in  the  case  of  OH  +  H2  and  OH  +  CH^  reactions,  there  is  a  sparsity  of 
kinetic  data  on  Reactions  (1)  and  (4)  over  the  relevant  combustion  and  plume 
temperatures.  Therefore,  k-j-k^  were  measured  in  the  temperature  range  of 
250-1000  K. 
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Experimental 


The  experimental  technique  of  flash  photolysis  resonance  fluorescence 
described  in  detail  in  Chapter  I  and  III  were  used  to  carry  out  these 


measurements. 


CO,  HC1 ,  CgHg,  and  CgHg  were  obtained  from  Matheson  Gas  Products.  They 
had  the  following  purity  levels:  CO,  99.99*;  HC1,  >99.99?* ;  CgH^,  >99.96%; 
and  CgHg,  >99.99%.  CO  was  passed  through  a  liquid  nitrogen  trap  before  use. 
HC1,  CgHg,  and  CgHg  were  subjected  to  degassing  cycles  before  use. 

Results  and  Discussion 

Experimental  measurements  of  k-j(T)  for  Reaction  (1)  are  listed  in 
Table  I  and  plotted  in  Figure  1.  As  seen  from  figure  1,  the  curvature  in 
the  Arrhenius  plot  is  quite  dramatic.  At  low  temperatures,  the  activation 
energy  for  this  reaction  is  near  zero;  at  higher  temperatures,  k^(T)  increases 
rapidly  with  temperature.  [It  is  worth  noting  here  that  visual  comparisons 
of  the  'degree'  of  curvature  in  Arrhenius  plots  for  several  reactions  can  be 
misleading;  for  a  reaction  where  the  low-temperature  rate  constant  itself 
depends  strong]  v  on  temperature,  e.g.,  OH  +  CH^  -*•  H^O  +  CH^,  substnatial  high- 
temperature  rate  constant  enhancement  may  be  perceived  as  only  moderate 
curvature.]  As  in  our  studies  of  OH  reactions  with  Hg  and  CH^,  our  measure¬ 
ments  of  k^(T)  fill  the  gap  between  500K  and  1000K  and  verify  the  adequacy 
of  the  empirical  forced  fit  expression  of  Baulch,  et  al.  Our  measurements 
of  kj(T)  above  900K  lie  slightly  above  the  solid  line  in  Figure  1;  it  would 
be  worthwhile  to  remeasure  and  extend  k^(T)  to  the  1200K  temperature  regime. 

Rate  coefficient  data  for  Reaction  (2)  are  listed  in  Table  II  and  plotted 
in  Figure  2.  Again,  the  fcnk2  vs.  1/T  plot  is  non-linear.  The  argeement 
between  previous  low- temperature  (T  <  298K)  measurements  and  present  results 
are  excellent.  However,  It  is  clear  that  extrapolation  of  low-temperature 
data  to  higher  temperatures  would  lead  to  large  errors. 
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Figure  1.  Plot  of  fcil^  vs. 

—  D.L.  Bauch,  D.O.  Drysdale,  J.  Duxbury,  and  S.  Grant, 

.  Evaluated  Kinetic  Data  for  High  Temperature  Reactions. 
Vol.  3,  Butterworths,  London,  p.  203-249. 
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Table  I.  Rate  Constant  Data  for  the 
Reaction  of  OH  with  CO. 


13 

k.  x  10  , 

Temperature  - 

(K)  cnrmolecule"  s~ 


298 

1.45  ± 

0.10 

511 

1.77  ± 

0.08 

607 

1.89  ± 

0.15 

703 

2.03  ± 

0.26 

800 

2.39  ± 

0.12 

922 

3.04  ± 

0.40 

992 


3.11  ±  0.40 


10  cm3 molecule-1  s 


1000 

TOO 

Figure  2.  Plot  of  itnk^  vs.  . 

I.W.M.  Smith  and  R.  Zellner,  JCS  Far.  Trans.  II,  7£,  1045  (1974). 

M.  Zahniser,  J.  Anderson,  and  F.  Kaufman,  Chem.  Phys.  Lett.  27_,  507  (1974). 

A.R.  Ravishankara,  G.  Smith,  Watson,  and  D.D.  Davis,  J.  Phys.  Chem. 

81,  2220  (1977). 
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Table  II..  Rate  Constant  Data  for  the 
Reaction  of  OH  with  HC1. 

k2  x  1012, 

Temperature  3  , 

(K)  cirrmolecule"  s~ 


252 

0.553 

±  0.036 

273 

0.656 

±  0.098 

298 

0.693 

±  0.051 

300 

0.645 

±  0.050 

328 

0.787 

±  0.083 

350 

0.850 

±  0.059 

357 

0.817 

±  0.054 

396 

0.930 

±  0.016 

402 

1.09 

±  0.10 

416 

1.21 

±  0.05 

438 

1.32 

±  0.07 

454 

1.15 

±  0.06 

500 

1.33 

±  0.16 

568 

1.56 

±  0.08 

576 

1.85 

±  0.12 

617 

1.88 

±  0.06 

656 

2.03 

±  0.22 

698 

2.58 

±  0.10 

752 

2.48 

±  0.22 

824 

2.69 

±  0.38 

932 

3.63 

±  0.46 

1000 

4.14 

±  0.37 
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Rate  constant  data  for  Reactions  (3)  and  (4)  are  tabulated  in  Table  III 
and  IV  and  plotted  in  Figures  3  and  4.  Prior  to  these  measurements ,  the  only 
reliable  rate  constants  existing  for  this  reaction  were  a  few  measurements 
at  and  slightly  above  room  temperature.  Once  again,  it  can  be  seen  from 
figures  3  and  4  that  straight-line  extrapolation  of  low-temperature  measure¬ 
ments  results  in  an  underestimation  of  the  bimolecular  rate  constants  at 
higher  temperatures;  concave  upward  curvature  is  agin  evident. 

The  chemical  kinetics  experiments  described  above  represent  valuable  new 
input  to  current  models  of  hydrocarbon  combustion  processes.  Reactions  (1) 
is  of  such  fundamental  importance  to  hydrocarbon  oxidation  chains  that 
accurate  characterization  of  these  k^(T)  functions  is,  in  itself,  a  goal 
of  considerable  significance.  Kinetic  estimates  for  Reactions  (1),  (3), 
and  (4)  pervade  hydrocarbon  combustion  modeling  efforts;  use  of  incorrect 
input  could  modify  the  bulk  observable  results  of  a  simulated  model  process 
to  an  extent  which  implies  that  the  assumed  model  structure  is  wrong— this 
may  be  the  case,  or  the  problem  could  be  in  the  inaccuracy  of  the  data  base. 

Yet  one  should  not  assume  that  the  value  of  these  results  (and  of  our 
developed  experimental  technique)  is  limited  to  the  numerical  achievements 
detailed  above.  Typically,  the  complexity  of  the  chemistry  involved  in  a 
combustion  and/or  plume  process  is  nearly  overwhelming.  In  recent  years 
it  has  become  increasingly  obvious  that  it  is  simply  not  possible  to  measure 
the  kinetic  parameters  for  all  contributing  elementary  reactions  in  a  com¬ 
plex  process.  Indeed,  if  such  a  data  base  did  exist,  computational  limita¬ 
tions  would  prevent  its  full  utilization.  Increasing  attention  is  thus 
being  focused  on  the  generation  of  'global'  kinetic  equations,  i.e.,  kinetic 
expressions  which  are  approximately  applicable  to  famililes  of  analogous 
chemical  reactions.  Such  a  viewpoint  then  argues  the  need  for  investigations 
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Table  III. 

Rate  Constant  Data  for  the 
Reaction  of  OH  with  ^Hg. 

Temperature 

00 

x  1012, 

cm  molecule  sec 

297 

0.259  ±  0.021 

400 

0.771  ±  0.076 

499 

1.58  ±  0.10 

609 

2.61  ±  0.33 

697 

3.65  ±  0.25 

800 

5.07  ±  0.34 
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Table  IV.  Rate  Constant  Data  for  the 
Reaction  of  OH  with  C^Hg. 

k.  x  1012, 

Temperature  ^ 

(K)  cm  molecule"  s~ 


298 

1.05 

+ 

0.04 

326 

1.48 

+ 

0.06 

378 

2.51 

+ 

0.20 

469 

3.37 

+ 

0.23 

554 

4.78 

+ 

0.34 

690 

8.78 

+ 

0.97 
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Figure  3.  Plot  of  Ink^  vs.  — y ■  . 
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which  emphasize  the  effects  of  chemical  structure  on  the  rates  of  reactivity; 
studies  aimed  at  generating  an  improved  understanding  of  the  general  tempera¬ 
ture  dependence  of  bimolecular  reaction  rates  must  be  pursued. 
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